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Abstract 

In [Ger05], the first author presented an efficient algorithm for computing involutive (and re- 
duced Grobner) bases. In this paper, we consider a modification of this algorithm which simpli- 
fies matters to understand it and to implement. We prove correctness and termination of the 
modified algorithm and also correctness of the used criteria. The proposed algorithm has been 
implemented in Maple. We present experimental comparison, via some examples, of performance 
of the modified algorithm with its original form described in [Ger05] and has been implemented 
in Maple too. In doing so, we have taken care to provide uniform implementation details for the 
both algorithms. 
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1. Introduction 



The most important algoritlimic object in computational algebraic geometry is Grohner 
basis. The notion of Grobner basis was introduced and an algorithm of its construction 
was designed in 1965 by Buchberger in his Ph.D. thesis [Buc65]. Later one [Buc79], he 
discovered two criteria for detecting some unnecessary, and thus useless, reductions that 
made the Grobner bases method a practical tool to solve a wide class of problems in poly- 
nomial ideals theory and in many other research areas of science and engineering [BW98] . 

In the 20s of the last century, French mathematician Janet [Jan20] developed a con- 
structive approach to analysis of certain systems of partial differential equations based 
on their completion to involution (cf. [SeilO]). The Janet approach was generalized to 
arbitrary polynomially nonlinear systems by American mathematician Thomas [Tho37]. 
Based on the related methods described in the book by Pommaret [Pom78], Zharkov 
and Blinkov [ZB96] introduced the notion of involutive bases in commutative algebra. 
The particular form of an involutive basis they used is nowadays called Pommaret basis 
[GB98,SeilO]. 

Gerdt and Blinkov [GB98] have proposed a more general concept of involutive bases 
for polynomial ideals and designed algorithmic methods for their construction. The un- 
derlying idea of the involutive approach is to translate the methods originating from 
Janet's results into polynomial ideals theory in order to construct involutive bases by 
combining algorithmic ideas of the theory of Grobner bases and those in the theory of 
involutive differential systems. In doing so, Gerdt and Blinkov [GB98] introduced the 
concept of involutive division. Moreover, they have applied the involutive form of Buch- 
berger's criteria. This led to a strong computational tool which is a serious alternative to 
the conventional Buchberger algorithm (note that any involutive basis is also a Grobner 
basis) . Apel and Hemmccke in [AH05] proposed two more criteria for detecting unneces- 
sary reductions in involutive basis computations (see also [Ger02]) which, in the aggregate 
with the criteria by Gerdt and Blinkov [GB98], are equivalent to Buchberger's criteria. 
The first author in [Gcr05] described an efficient algorithm to compute invohitive and 
Grobner bases using all these criteria. We refer to Seller's book [SeilO] for a comprehen- 
sive study and application of involution to commutative algebra and geometric theory of 
partial differential equations. 

In this paper, we propose a variant of Gerdt's algorithm [Ger05] for constructing invo- 
lutive bases. This algorithm seems to be simpler both to understand and to implement. 
We prove correctness for the new version of algorithm which includes the criteria and its 
termination. We have implemented this algorithm in Maple and we compare here its per- 
formance with our implementation of the original Gerdt's algorithm via some examples. 

The structure of the paper is as follows. Section 2 contains the basic definitions and 
notations related to theory of involutive bases, and a short description of the algorithm 
for their computing proposed in its initial form [GB98]. In Section 3, we briefly present 
Gerdt's algorithm from [Ger05] to compute involutive bases. Section 4 is devoted to 
description of the modified algorithm together with a proof of its correction and termi- 
nation. In Section 5, we compare, by using some benchmarking examples, performance 
of the last algorithm with our implementation of Gerdt's algorithm done also in Maple. 
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2. Preliminaries 



In this section, wc recall some; basic definitions and notations from the theory of 
involutive bases which are used in the paper, and briefly describe the initial algorithm 
proposed in [GB98] for computing such bases. 

Let K he a field and R = K[xi, . . . , x„] the polynomial ring in the variables Xi, . . . ,Xn 
over K. Below, we denote a monomial x"^ • • • e i? by x" where a = (ai , . . . , a„) e N" 
is a sequence of non-negative integers. A monomial ordering on R used in theory of 
Grobner and involutive bases is a total order -< on the set of all monomials satisfying the 
following two properties: 

• It is multiplicative; i.e., x" -< x^ implies x"+''' -< x^'^'^ for all q;,/3,7 G N", 

• The constant monomial is the smallest; i.e., 1 ^ x" for all a € N". 

A typical example of a such monomial ordering is the pure lexicographical ordering, 
denoted by -<iex- For two monomials x",x^ G i? x" -<;ex if the left most nonzero 
entry of /? — a is positive. 

Let / = (/i, . . . , /fc) be the ideal of R generated by the polynomials fi,. . . ,fk G R. 
Also let / e i? and ^ be a monomial ordering on R. The leading monomial of / is the 
greatest monomial (with respect to -<) appearing in /, and wc denote it by LM(/). If 
F C i? is a set of polynomials, we denote by LM{F) the set {LM(/) | f G F}. The 
leading coefficient of /, denoted by LC(/), is the coefficient of LM(/). The leading term 
of / is LT(/) = LC(/)LM(/). The leading term ideal of / is defined to be 

LT(/) = (LT(/) I /e/). 

A finite set G = {gi, . . . , gk} C / is called a Grobner basis of / w.r.t. -< if LT(7) = 
(LT(5i), . . . ,LT(£ffe)). For more details, we refer to [BW93], pages 213-214. 

We recall below the definition of a special kind of involutive bases, namely, Janet basis. 
For this purpose, we describe first the notion of an involutive division [GB98] which is a 
restricted monomial division [Ger05] that, together with a monomial ordering, determines 
the properties and the structure of an involutive basis. This makes the main difference 
between involutive bases and Grobner bases. The idea is to partition the variables into 
two subsets of multiplicative and non-multiplicative variables, and only the multiplicative 
variables can be used in the divisibility relation. 

Definition 1. ([GB98]) An involutive division £ on the set of monomials of R is given, 
if for any finite set U of monomials and any u €U, the set of variables is partitioned into 
the subsets of multiplicative Mc{u, U) and non-multiplicative NMc{u, U) variables, and 
the following three conditions hold: 

• u,v gU, ujC{u, U) n vC{v, U) T^^^^ue vjC{v, U) or v G ujC{u, U) 

• V gU, V e uC{u, U) =^ C{v, U) C C{u, U) 

• uGV imAV ^ C{u, U) C £(u, V) 

where C{u,U) denotes the monoid generated by the variables in Mc{u,U). If v G 

uC{u,U) then u is called C— (involutive) divisor of v and the involutive divisibility re- 
lation is denoted by u\cv. If v has no involutive divisors in a set U, then it is called 
£— irreducible modulo U. 

There are involutive divisions based on the classical partitions of variables suggested 
by Janet [Jan20] and Thomas [Tho37]. In this paper, we are interested only in Janet 
division [GB98]. 
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Example 2. (Janet division) Let U he a. finite set of monomials and u £ U . The 

variable xi is multiplicative for u, if u has the maximum degree of xi in U. For z > 1, 
the variable Xi is multiplicative for u if it has the maximum degree of Xi in the set 
{v ^ U \ Aegj{v) = Aegj{u),l < j < i}, where degj{v) denotes the degree of Xj in a 
monomial v. 

Proposition 3. ([GB98]) Janet division satisfies the conditions in Definition 1 and is 
noetherian. 

Throughout this paper £ is assumed to be the Janet division. Now, we define an invo- 
lutive basis. 

Definition 4. ([Ger05]) Let J C i? be an ideal and -< a monomial ordering on B. Let 
£ be an involutive division. A finite set G C / is an (£-)involutive basis for / if for all 
/ e / there exists g €G such that LM{g)\cLM{f). 

From this definition and from that for Grobner basis [Buc65,BW93] it follows that an 
involutive basis for an ideal is its Grobner basis, but the converse is not always true. 

Remark 5. By using an involutive division in the conventional division algorithm for 
polynomial rings, we obtain an involutive division algorithm. If F is a finite polynomial 
set and C is an involutive division, then we use NF£(/, F) to denote the remainder of / 
on involutive division by F. 

The following theorem provides an algorithmic characterization of involutive basis 
for a given ideal which is an involutive analogue of the Buchberger characterization of 
Grobner basis. 

Theorem 6. ([GB98]) Let I C R be an ideal, -< a monomial ordering on R, and C an 
involutive division. Then a finite subset G G I is an involutive basis for I if for each 
f eG and each x € NMc(LM{f),LM{G)), we have NFc{xf,G) = 0. 

Based on this theorem, one can design an algorithm to compute involutive bases. We 
recall here the InvBasis algorithm from [GB98] (see also [Ger05]) which computes a 
minimal involutive basis for an ideal. An involutive basis G is called minimal if for any 
other involutive basis G the inclusion LM(G) C LM(G') holds. A minimal involutive basis 
exists and being monic and involutively autoreduced, i.e. satisfying 

(VfleG) [g = NFc{g,G\{g})] 

is unique for a given ideal, a monomial ordering and a constructive involutive division 
(see [GB98] for the definition of constructivity) . 



4 



Algorithm 1 InvBasis 

Input: F, a set of polynomials; £, an involutive division;, -<, a monomial ordering. 
Output: a minimal involutive basis for (F) w.r.t £ and -< 
Select f e F with no proper divisor of LM(/) in LM{F) 

G := {/}; 
Q:=F\G; 
while Q 7^ do 

Select and remove p £ Q with no proper divisor of LM{p) in LM{Q); 
h:= NFc{p,G); 
if ft. 7^ then 

for g E G which LM{g) is properly divisible by LM(/i) do 

g :-gu{5}; 

G:=G\{g}; 
od 

G:=GU{/i}; 

Q:=QU{xg\gGG,xe NMc{LM{g), LM(G))}; 

fi 
od 

Return (G) 



Here, as it has been mentioned in [Ger05], in comparison to the algorithm in [GB98] 
another selection strategy is used. By this strategy, a polynomial whose leading monomial 
has no proper divisor is chosen. However, this algorithm is not cfScicnt in practice, it 
processes the repeated prolongations and does not use any criterion to avoid the unnec- 
essary reductions. For this purpose, we consider in the next section, the Gerdt' algorithm 
from [Ger05] as an improvement of InvBasis. 

3. Gerdt's algorithm 

The Gerdt's algorithm [Ger05] improves the above algorithm InvBasis. It uses the 
involutive form of Buchberger's criteria to avoid unnecessary reductions and avoids the 
repeated processing of non-multiplicative prolongations. In order to explain the structure 
of this algorithm, we recall some more definitions and notations. 

Definition 7. ([Gcr05]) An ancestor oi a polynomial f <E F C i?\{0}, denoted by anc(/), 
is a polynomial g £ F oi the smallest deg(LM(5')) among those satisfying LM(/) = 
uLM{g) where u is either the unit monomial or a power product of non-multiplicative 
variables for LM(fif) and such that NF£(/ -ug,F\ {/}) = if / 7^ ug. 

This additional information for an clement in a polynomial set is useful to avoid un- 
necessary reductions specially by applying the adapted Buchberger's criteria (see below). 

Proposition 8. ([AH05]) Let I C R be an ideal and G C I be a finite set. Let also -< 
be a monomial ordering on R and £ an involutive division. Then G is an C— involutive 
basis for I if for all f & G one of the following conditions holds (we use "C " to denote 
proper conventional division): 
(1) for all X G NMc(LM{f),LM{G)) the equality mc{xf,G) = holds. 
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(2) There exists g & G with LM(p)|£LM(/) satisfying one of the following: 
(Ci) LM(anc(/))LM(anc(.g)) = LM(/) 
(C2) lcm(LM(anc(/)),LM(anc(5))) □ LM(/) 
(C3) There exists t € G such that: 

• lcm(LM(t),LM(/)) C lcm(LM(anc(/)), LM(anc(g))) 

• lcm(LM(i),LM(5)) □ lcm(LM(anc(/)), LM(anc(5))) 

(C4) There exists t G G computed before anc(/) and y e N M c{lM.{t) ,IM.{G)) s.t. 

• ylM{t) = LM(/) 

• lcni(LM(anc(/)),LM(anc(t))) CLM(/). 

Now we can present Gerdt's algorithm for computing involutive bases in which we 
associate to each polynomial /, the triple p = {f^g^V} where / = poly(p) is the poly- 
nomial itself, g — anc(p) is its ancestor and V = NM{p) is the list of non-m\iltiplicativc 
variables of / which have been already processed in the algorithm. If P is a set of triples, 
we denote by poly(P) the set {poly(p) | p G P}. If no confusion arises, we may refer to 
a triple p instead of poly(p), and vice versa. 



Algorithm 2 Gerdt 

Input: F, a set of polynomials; £, an involutive division; -<, a monomial ordering 

Output: a minimal involutive basis for {F) w.r.t C and ^ 
Select f £ F with no proper divisor of LM(/) in LM(P) 

T:={{f,fM}; 
Q:={{g,g,0} \qeF\{f}}- 
Q :=HeadReduce((3, T, £, -<); 
while Q ^ do 

Select and remove p ^ Q with no proper divisor of LM(poly(p)) in LM(poly(Q)); 
if poly(p) = anc(p) then 

for g e T with LM(poly(p)) C LM(poly(g)) do 
Q:=QU{q}; 
T:=T\{q}; 
od 

fl 

h :=TailNormalForm(p, T, jC, 
T:=TU{{h,anc{p),NMip)}}; 

for g e T and x e 7VAf£(LM(poly(g)), LM(poly(T))) \ NM{q) do 
Q:=QU {{x poly(g), anc(g), 0}}; 

NM{q) := NM{q) D TVM^ (LM (poly (g) ), LM (poly (T))) U {x}; 
od 

Q :=HeadReduce((5, T, £, ^); 
od 

Return (poly(T)) 



This algorithm includes three subalgorithms Head Reduce, HeadNormalForm and 
TailNormalForm which we present below. 
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Algorithm 3 HeadReduce 



Input: Q and T, sets of triples; £, an involutive division; -<, a monomial ordering 
Output: The Q of triples whose polynomials are £-head reduced modulo T 

S:=Q 
Q:= 0; 

while Sjt<D do 

Select and remove p G 5 

h :=HeadNormalForm(p, T, £); 

if hy^O then 

if LM(poly(p)) ^ LM{h) then 

Q ■.= QU{{h,h,9}}; 
else 

Q:=QU {p}; 

fi 

else 

if LM(poly(p)) = LM(anc(p)) then 
S:=S\{q€S\ anc(g) = poly(p)}; 

fi 

fl 
od 

Return (Q) 



In the below subalgorithm HeadNormalForm the Boolean expression Criteria(p, g) 
is true if at least one of the four criteria in Proposition 8 holds for p and g. 

Algorithm 4 HeadNormalForm 

Input: T, a set of triples; p, a triple; £, an involutive division; a monomial ordering 
Output: £-head normal form of poly(p) modulo T 

h := poly(p); 
G := poly(r); 

if LM(/i) is £-irreducible modulo G then 

Return {h) 
else 

Select .9 e G with LM(poly(£f))|£LM(/i); 
if LM{h) ^ LM(anc(p)) then 
if Criteria(p, 5) then 
Return (0) 

fi 
else 

while hj^O and LM{h) is £-reducible modulo G do 
Select g€G with LM{g)\cLM{h)] 



od 

fi 

fi 

Return (h) 
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Algorithm 5 TailNormalForm 

Input: T, a set of triples; p, a triple, which is £-head reduced modulo T; £, an involutive 

division; -<, a monomial ordering 
Output: £-Normal form of poly(p) modulo T 

h := poly(p); 
G:=poly(r); 

while h has a term t which is /^—reducible modulo G do 

Select 5 e G with LM{g)\ct; 

h:=h- gjj^^; 
od 

Return (h) 



4. Modified algorithm 

In this section, we present the following variant of Gerdt's algorithm for computing 
involutive bases. 

Algorithm 6 VarGerdt 

Input: F, a set of polynomials; £, an involutive division; -<, a monomial ordering 
Output: a minimal £— involutive basis for (F) 

Select f G F with no proper divisor of LM(/) in LM{F) 

T :={{/,/, 0}}; 
Q:={{q,q,n \q€F\{f}}; 
while Q ^ do 

Select and remove p & Q with no proper divisor of LM(poly(p)) in LM(poly((5)); 

h :=NormalForm(p, T, C, -<); 

if /i = and LM(poly(p)) = LM(anc(p)) then 

Q:={q&Q\ anc(g) ^ poly(p)}; 
else 

if LM(poly(p)) ^ LM(/i) then 

for g G T with LM(poly(/i)) C LM(poly(g)) do 

Q :=QU{q}; 

T:=T\{q}; 
od 

T:=TU{{h,h,9}}; 
else 

T:=TU{{h,ancip),NM{p)}}; 

fi 

for g e T and x G NMc{hM{poly{q)), LM(poly(T))) \ NM{q) do 
Q:=QU{{x poly(<7),anc(<7),0}}; 

NM{q) := iVM(g) n iVM£(LM(poly(g)), LM(poly(T))) U {x}; 
od 

fi 
od 

Return (poly(T)) 
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This variant seems to be simpler both to understand and to implement. We also present 
here a proof of its correctness and termination including the correctness proof for the 

used criteria. 

The subalgorithm NormalForm returns the £— normal form of a polynomial w.r.t. 
a given set of polynomials (sec below). Moreover, this subalgorithm detects some unnec- 
essary reductions using the involutive form of Buchberger's criteria. It is worth noting 
that the first author and Yanovich [GY06] have shown that the use of C3 and C4 criteria 
(see Proposition 8) may notably slow down the computation of involutive bases. That 
is why, we use only Ci and C2 criteria. Below, we give a simple proof of correctness for 
these criteria. 

Lemma 9. Let I (Z R be an ideal, -< a monomial ordering on R and C an involutive 
division. Let T c I be the last computed basis during a computation of an involutive 
basis for I, and f e I. Then NFc{f,T) = if there exists q gT with LM{q)\c'LM{f) 
satisfying one of the following conditions: 

(Ci) LM(anc(/))LM(anc(9)) = LM(/) 

(C2) lcm(LM(anc(/)),LM(anc(9))) C LM(/). 

Proof. Suppose that (/, q) satisfies Ci. Then, for lcm(LM(anc(/)), LM(anc((7))) two cases 
are possible: If it is a proper divisor of LM(anc(/))LM(anc(g)) = LM(/), there exists a 
monomial s ^ 1 such that LM(/) = slcm(LM(anc(/)), LM(anc(g))). Let f = u - anc(/), 
q = V- anc(g) and LT(/) = t ■ LT(g) for some monomials u and v and some term t. Thus, 
we can write f — tq = u - anc(/) — tv ■ a,nc{q) = s{u' ■ anc(/) — v' ■ anc(g)) where u = su' 
and vt = sv' for some monomials u',v'. Since, lcm(LM(anc(/)), LM(anc(g))) is a proper 
divisor of LM(/), then H'-anc(/)— ti'-anc((7) has been computed before / (see the selection 
strategy for choosing polynomials in the algorithm), and therefore, it has a standard 
representation w.r.t. T. This implies that f ~tq has also a standard representation w.r.t. 
T, and NF£(/, T) = 0. As the second case, if lcm(LM(anc(/)), LM(anc(g))) is equal to 
LM(anc(/))LM(anc(g)), then by the Buchberger's first criterion, we can conclude that 
u' ■anc{f)—v' ■anc{q) has a standard representation w.r.t. T, and this proves the assertion 
like the first case. 

Now, if (/, q) satisfies C2, the proof is similar to that for the above first case. □ 

,Just as above, for two polynomials p and q, the Boolean expression Criteria(p, q) is 
true if they satisfy at least one of the criteria Ci or C2, and false otherwise. 
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Algorithm 7 NormalForm 

Input: p, a triple; T, a set of triples; £, an involutive division; a monomial ordering 
Output: £-Normal form of p modulo T 

h := poly(p); 
G := poly(r); 

while h has a term t which is £— reducible modulo G do 
Select g eG with LM{g)\ct; 
if t = LT(poly(p)) and Criteria(/i., 3) then 

Return (0) 
else 

h:= h- Qij^^; 

R 
od 

Return (h) 



Theorem 10. The VarGerdt algorithm computes a minimal involutive basis for the 
ideal generated by the input polynomial set. 

Proof. Lot C -R be a finite set of polynomials taken as the input of VarGerdt 
algorithm. Let ^ be a monomial ordering on R and £ bo a constructive nocthcrian 
involutive division [GB98]. Let G be a a polynomial set that is the output of the al- 
gorithm for F. Since (by the structure of the algorithm) for each f € G and each 
X € NMc(L'M{f), LM(G)), we examine xf (we either compute its £— normal form w.r.t. 
G or we discard it by the criteria), it follows that NFcixf.G) = 0. Therefore, G is an 
involutive basis for (F) by Theorem 6. Moreovc;r, this basis is minimal by the definition 
of this concept due to the structure of the VarGerdt algorithm and constructivity of 
C. 

The termination of VarGerdt algorithm follows from the noetherianity of C, see 
Proposition 3. □ 

Remark 11. The algorithm VarGerdt as well as Gcrdt's algorithm not necessarily 
outputs an involutively tail autoreduced basis. If the output is G, then, in accordance to 
the definition that follows Theorem 6 of Section 2, the tail autoreduction is provided by 
the command 

(V<?€G) [g:=NFc{g,G\{g})]. 

5. Experiments and results 

We have implemented both algorithms VarGerdt and Gerdt in Maple 14 ^ . For an 
efficient implementation of Gerdt algorithm, we refer to [BGC03]. It is worth noting 
that, in this paper, we are willing to compare the structure and behavior of VarG- 
erdt and Gerdt algorithms on the same platform. Therefore, we do not compare our 
implementations with [BGC03]. 

To compare the behavior of these algorithms, we used some well-known examples 
from the the collection [EM] have been widely used for verification and comparison of 



The Maple code of our programs and examples is available on the Web page http://invo.jinr.ru/ 
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different software packages for construction of Grobner bases. In our benchmarking with 
these examples we apphed the following selection strategy for the polynomials in Q. An 
element at the initialization step to be inserted in set Q as well as that selected from 
Q in the main loop has the lowest leading monomial with respect to the degrce-reverse- 
lexicographic ordering. In the case of several such polynomials that is selected whose 
ancestor has been examined earlier than the ancestors of the other polynomials with the 
identical leading monomials. 

The results are shown in the following tables where computation was performed on 
a personal computer with 2.33 GHz, 2xIntel(R) Xeon(R) Quad core, 16 GB RAM and 
64 bites under the Linux operating system. The computation was done over Q and the 
monomial ordering is always the degrec-rcversc-lexicographical one. 
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655916.5 


476 


152 


18 





10 



The time (resp. memo., reds., Ci and C2) column shows the consumed CPU time in 
seconds (resp. amount of megabytes of memory used, number of reductions to zero, the 
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number of polynomials removed by Ci and C2 criteria) by the corresponding algorithm. 
The seventh column indicates the number of polynomials eliminated by Rewritten crite- 
rion: in VarGerdt algorithm, as well as in Gerdt algorithm, when the polynomial part 
of a triple {/, g, V} reduces to zero, and f = g, then we can discard any polynomial whose 
ancestor is /. We call this criterion the Rewritten criterion. The last column shows the 
largest degree of intermediate polynomials treated during the computation of involutive 
bases. 

A comparison of the timing columns in the above tables and our test for some other 
examples shows that VarGerdt is more efficient and stable than Gerdt. Indeed, the 
main difference between these algorithms is that in Gerdt algorithm, we do first the 
head reduction of all non-multiplicative prolongations (which have not been examined 
yet) and then we choose one of them to complete its involutive reduction to the full 
normal form, while in VarGerdt we choose a polynomial (among the non-multiplicative 
prolongations), and we compute its full normal form. 

In Gerdt 's algorithm the head involutive reduction of the whole set of non- multiplicative 
prolongation is done in order to provide a platform for matching a good selection strat- 
egy. We refer to paper [GB07] for details on heuristically good selection strategies for 
the Janet division. The choice of such selection strategy as well as the use of proper data 
structures for fast search of involutive divisor [Gcr05] plays a key role in providing a high 
computational cifficicincy of the involutive algorithms. In this paper we do not consider 
these important aspects of the practical construction of involutive bases and compare 
two algorithms experimentally for their simplest and identical for the both algorithms 
implementation. 
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